Observing Trans- Planckian Signatures in the Cosmic Microwave Background 
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We examine the constraints cosmological observations can place on any trans-Planckian correc- 
tions to the primordial spectrum of perturbations underlying the anisotropics in the Cosmic Mi- 
crowave Background. We focus on models of trans-Planckian physics which lead to a modulated 
primordial spectrum. Rather than looking at a generic modulated spectrum, our calculations are 
based on a specific model, and are intended as a case study for the sort of constraints one could hope 
to apply on a well-motivated model of trans-Planckian physics. We present results for three different 
approaches - a grid search in a subset of the overall parameter space, a Fisher matrix estimate of 
the likely error ellipses, and a Monte Carlo Markov Chain fit to a simulated CMB sky. As was seen 
in previous analyses, the likelihood space has multiple peaks, and we show that their distribution 
can be reproduced via a simple semi- analytic argument. All three methods lead to broadly similar 
results. We vary 10 cosmological parameters (including two related to the trans-Planckian terms) 
and show that the amplitude of the tensor perturbations is directly correlated with the detectability 
of any trans-Planckian modulation. We argue that this is likely to be true for any trans-Planckian 
modulation in the paradigm of slow-roll inflation. For the specific case we consider, we conclude 
that if the tensor to scalar ratio, r 0.15, the ratio between the inflationary Hubble scale if, and 
the scale of new physics M has to be on the order of 0.004 if the modulation is detectable at the 2cr 
level. For a lower value of r, the bound on H/M becomes looser. 



I. INTRODUCTION 

Conventional cosmology is based upon a four dimen- 
sional spacetime whose evolution is governed by the Ein- 
stein field equations. This picture certainly breaks down 
at the Planck scale due to quantum corrections to the 
gravitational action. Furthermore, stringy effects may 
change our understanding of spacetime at energies one 
or two orders of magnitude below the Planck scale. Un- 
fortunately, both the Planck and string energy scales are 
entirely inaccessible to direct experimentation. However, 
thanks to the dramatic expansion of the universe during 
the inflationary era, it is possible that cosmology pro- 
vides an unexpected window into physics at the shortest 
length scales jlj-jl6j. 

If we assume the universe has always been matter dom- 
inated and is now around 13 billion years old, the current 
horizon size is roughly 40 billion lightyears, or 2 x 10^^ 
Planck lengths. If inflation ended at the GUT scale 
(~ IQ-'^^GeV) with efficient reheating, the universe has 
grown roughly 4 x 10^^ times larger since reheating, and 
by at least a similar amount during inflation, in order to 
ensure that the universe is homogeneous on scales equal 
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to the present horizon size. Consequently, all distances 
significantly smaller than galactic scales in the current 
universe were necessarily sub-Planckian at some point 
during inflation. If the inflationary period continues for 
longer than the minimum time necessary to solve the 
cosmological flatness and horizon problems, all presently 
visible scales must have crossed the Planckian boundary 
during inflation. Conversely, without inflation all astro- 
physical scales are mapped to lengths larger than the 
Planck length at all times since the age of the universe 
was of order the Planck time. 

On the face of it, the fact that astrophysical scales in 
the present universe correspond to sub-Planckian scales 
during the inflationary epoch is simply a mathemati- 
cal curiosity. However, inflation generates perturbations 
through the exponential stretching of quantum fluctua- 
tions. The textbook calculation of the primordial per- 
turbation spectrum implicitly ignores the consequences 
of the breakdown of spacetime at some new fundamen- 
tal lengthscale. For a clear exposition of the standard 
treatment of perturbations, including a discussion of the 
effect of boundary conditions, see 17]. In string theory 
the effective minimum length can be significantly larger 
than the Planck length, by perhaps one or two orders of 
magnitude 18, 19, 20]. In this case, the ratio between 
the inflationary Hubble scale, and the scale of new 
stringy physics, M can be as large as 10~^. 
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Dimensional analysis alone suggests that any modifi- 
cations to the perturbation spectrum induced by includ- 
ing this minimum length will be a function of (i^/M)^, 
where p is, of course, unknown. If p is unity, then there 
is at least a corner of parameter space where we might 
expect potentially observable modifications to the spec- 
trum. However, for larger values of any modifications 
to the spectrum are most likely unobservable. String 
theory and other candidate unified models are not suf- 
ficiently mature to permit an ab initio analysis of per- 
turbation generation, and predict the value of p. Con- 
sequently, all studies of trans-Planckian modifications 
to the infiationary perturbation spectrum require some 
ansatz which, while motivated by our understanding of 
stringy / Planckian physics, introduces a degree of ar- 
bitrariness into the calculation. This ansatz necessarily 
modifies either the initial conditions or the evolution of 
the perturbations. Generically, modifications to the ini- 
tial conditions produce changes in the spectrum on the 
order of H/M,^ whereas modifications to the evolution 
result in changes on the order of {H/M)'^. It is worth 
noting, however, that no calculation has predicted that 
the modifications vanish completely. 

Normally, this lack of theoretical consensus would de- 
ter detailed studies of the observational consequences of 
any given model of trans-Planckian physics. However, 
even a small possibility that the fingerprints of Planck- 
ian or stringy physics can be found on the sky provides 
more than enough motivation for us to take a careful 
look at what constraints can be established with present 
data or future measurements. In this paper, we focus on 
one model of trans-Planckian physics, first suggested by 
Danielsson in the case of de Sitter space 10] and gener- 
alized to non-de Sitter backgrounds by Easther, Greene, 
Kinney and Shiu 11]. There is nothing intrinsically spe- 
cial about this model - our analysis is effectively a case 
study of what might be achieved when one has a specific 
model to test. In particular, we do not enter into the 
debate as to whether or not such "truncated a-vacua" 
are self-consistent with respect to being well-defined in 
the ultraviolet limit or with respect to their suitability 
as initial conditions for infiation. We simply consider the 
model as a reasonable case study for investigating the 
discriminatory power of forthcoming CMB observations. 

The one general statement we can make is that it is far 
more likely that trans-Planckian effects will have an ob- 
servable impact on the spectrum if infiation takes place 
at a relatively high energy scale. In principle, infiation 
can occur at almost any scale high enough to allow for 
baryogenesis and nucleosynthesis in the post-infiationary 
universe - and this can be far below the GUT scale. How- 
ever, since we only expect a detectable signal if H/M 
is comparatively large, we are pushed into the corner 
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of parameter space where M is small (compared to the 
Planck scale) while H is roughly GUT-scale. A small 
M (two orders of magnitude below the Planck scale, 
say) is entirely consistent with string theory, whereas 
H is effectively a free parameter. An upper bound on 
H arises from the limits on the tensor component to 
the primordial perturbation spectrum - which is poten- 
tially observable via both the temperature anisotropics 
in the CMB, and the B-mode of the CMB polarization. 
(A rough upper limit on H from the WMAP data is 
H < 6 X 10^^ GeV.) It is very easy for infiation to hap- 
pen at a low enough scale to ensure that a primordial 
tensor spectrum is effectively unobservable by any con- 
ceivable experiment [22, IIM ISJ? l23 • However, if we were 
to observe a contribution to the CMB that appears to be 
related to Planckian or stringy physics, there is a strong 
likelihood that we would be in the region of parameter 
space where H was large enough in order to ensure that 
the tensor contribution was also detectable. 

There are important consequences to this qualitative 
correlation between trans-Planckian contributions to the 
CMB and the tensor perturbations. First, it further 
raises the stakes for high precision measurements of the 
CMB. We already know that an observable tensor spec- 
trum fixes the infiationary energy scale and rules out 
some of infiation's competitors, such as the ekpyrotic 
and pre-Big Bang scenarios |23, |2^ , and that single field 
models of infiation must satisfy a set of consistency condi- 
tions 29]. However, on top of these achievements, deter- 
mining H would let us put direct constraints on the scale 
of new physics in any specific model of trans-Planckian 
physics, since M would be the only remaining free pa- 
rameter. If H turns out to be at the high end of the 
permitted range, it increases the likelihood that we will 
be able to see any trans-Planckian signal that is there 
to be found. Finally, one us [RE] recently argued that a 
high value of H during infiation is naturally correlated 
with a scalar spectrum with a significant running index 
{dns/d\nk ^ 0) [s^. We do not consider models with 
substantially broken scale invariance in this paper, but do 
incorporate a non-zero running into our analysis below. 

From a practical perspective, if we perform a fit to 
CMB data with a model that includes trans-Planckian 
effects, it is important to include the tensor modes in 
any analysis, since the same parameter {H) that governs 
their amplitude and scale dependence is also likely to 
appear in the trans-Planckian corrections. Conversely, 
simply observing an oscillatory spectrum would not be 
enough to prove that Planck / string scale physics had 
left an imprint on the CMB - one can produce the same 
result purely with physics at lower energy |3i|, |^ . 

This problem has been addressed in several previous 
analyses. The first computation of trans-Planckian mod- 
ifications to the spectrum was |^, which used a Fisher 
matrix calculation to predict that the Planck mission will 
put stringent bounds on a oscillatory component in the 
fundamental spectrum. However, this analysis marginal- 
ized over a small subset of cosmological parameters, and 
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the use of a Fisher matrix in this context is problematic 
|33 |. Subsequent calculations have been less optimistic. 
In particular, argued that any reasonable modulation 
to the spectrum is most likely unobservable, although we 
will argue that their pessimism is probably misplaced. 
Okamoto and Lim |23 looked at the constraints on a 
general superimposed oscillation to the spectrum. They 
varied several cosmological parameters as well as those 
describing the oscillation, and concluded that values of 
H/M>0.005 might be reached with a cosmic variance 
limited survey, but ignored the role of tensor modes. 
Martin and Ringeval, |3^ |37|, |3S] have studied models 
with independent parameters governing the amplitude 
and frequency of the trans-Planckian modulation. Their 
analysis is distinguished by the possibility that very high- 
frequency oscillations with non-negligible amplitude are 
considered, with the frequency of the oscillations produc- 
ing the primary constraint on H/M. In such a case, they 
conclude that very strong constraints are possible on the 
value of H/M^ even when restricted to current data. 

This paper is organized as follows: In Sec. we dis- 
cuss the specific model of trans-Planckian physics we con- 
sider in this paper. The modulation to the power spec- 
trum is correlated to the tensor/scalar ratio, which we 
argue is well motivated because of its relationship to a 
/^-independent boundary condition. In Sec. IIIIl we dis- 
cuss two complementary surveys of the trans-Planckian 
parameter space: First, we consider a grid-based search 
of the parameter space, which is computationally ineffi- 
cient but enables a broad sweep of the space. Second, we 
consider a Fisher-matrix based analysis, which handles 
large numbers of parameters in an efficient way, but does 
not sample the global properties of the likelihood surface. 
In Sec. IIVI we describe a Markov Chain Monte Carlo 
analysis, which constitutes a complementary method for 
exploring the parameter space. All three analyses are 
in good agreement. Comments and conclusions are pre- 
sented in Sec. 



II. THE MODEL 

Quantum fiuctuations during infiation produce two 
fiuctuation spectra: the curvature perturbation in the 
comoving gauge 7^, and the two polarization states of 
the primordial tensor perturbation, /i+ and h^. We pa- 
rameterize these power spectra by 
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where A^(/c^) is a normalization constant, and is some 
pivot wavenumber. The running, dn/dhik^ is defined by 
the second derivative of the power spectrum, dn/dhik = 
d^A'^/dlnk^, for both the scalar and the tensor modes. 



This parameterization gives the definition of the spectral 
index. 
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for the scalar modes, and 
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for the tensor modes. In addition, we re-parameterize the 
tensor power spectrum amplitude, A^(/c^), by the "ten- 
sor/scalar ratio r", the relative amplitude of the tensor- 
to-scalar modes, given by 



A|(fc.) 



(5) 



For a single slowly rolling infiaton, nt = — r/8, so we can 
reduce the number of parameters by fixing nt once r is 
chosen. We choose the pivot-scale k^ = 0.002 Mpc~^ at 
which to evaluate ng, r and A, where A{k^) and A'^{k^) 
are related [see 39j through 



A^(fc^) 2.95 X 10-^ A{k^). 



(6) 



Note that in this definition, the trans-Planckian parame- 
ter e, which appears in Eq. is simply the first slow-roll 
parameter [see e.g. il7j . 
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where Mpi = (STrG)-^/^ ^ mpi/V^ = 2.4 x 10^^ GeV 
is the reduced Planck energy; it is related to the tensor- 
to-scalar ratio by r = 16e. The definition in terms of 
i^(^), called the Hubble slow roll expansion, is exact, 
and the expression in terms of the infiationary potential 
V(^), called the potential slow roll expansion, is an ap- 
proximation valid in the limit of a slowly rolling field, 
<C V {(/)). The primordial power spectra A'^ and A^ 
are the underlying spectra which are modulated by ef- 
fects from Planck-scale physics, described below. 

The specific model of trans-Planckian physics we are 
working with 13, Til assumes that rather than reducing; 
to the Minkowski mode function Uk in the infinite past, 
the perturbations are matched at some finite time, cor- 
responding to the moment (different for each comoving 
perturbation) when their physical wavelength exceeds the 
minimum length, 1/M. Following pjj, we can compute 
the modified spectrum; 



BD 



(8) 



where BD refers to the spectrum obtained with the 
conventional Bunch-Davies vacuum. In this case the 
scalar and tensor spectra are both modified by the same 
amount. The C± are 11] 

^ I f iyc 

^- = -o exp 

2 VI - \yc 



- ]uk{yc)- (10) 
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Here is a rescaled "time", defined hy y = k/{aH) |40j . 
so that in the absence of a new scale y decreases from 
+00 to 0, with horizon crossing taking place (by defini- 
tion) when y = 1. In this language, ydk) is the criti- 
cal time for the k-th mode, at which its physical scale 
becomes larger than the minimum length. This is a k- 
dependent quantity, since in all slow roll models, H is 
a (slowly varying) function of time, and in general |Tl| . 
yc cx k^ where e is the usual slow roll parameter. This 
relationship is exact when e is expressed in the Hubble 
slow roll expansion, and approximate when the potential 
slow roll form is used. Consequently, 



(11) 



where M is the scale of new physics, and the starred 
quantities again refer to some fiducial value of k. The 
quantity Uk{yc) is the value of the quantum mode func- 



tion (given in the short-wavelength limit by Uk 
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evaluated when the wavelength of the mode is equal to 
a fixed physical cutoff length k/a = M, or equivalent ly 
y = yc(k). See Ref. Tl] for a detailed discussion. 

Keeping only the lowest order term in we can 

write the modified spectrum in a more transparent form. 



^^l . ( 2y, 
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where is a phase factor quantifying our lack of a priori 
knowledge about the physical "pivot-scale" k^. 

We note that the functional form of the modification 
above is determined not by the details of the boundary 
condition, but by the behavior of the background, in par- 
ticular the changing ratio of the infiationary horizon size 
to the (fixed) fundamental scale. Such a sinusoidal modu- 
lation likely to be a generic feature of any trans-Planckian 
ansatz that changes the initial conditions for a given per- 
turbation, so long as the boundary condition itself has no 
intrinsic dependence on /c.^ Since this effect depends on 
which gradually decreases during the course of the in- 
fiationary epoch, the original spectrum is multiplied by 
a (small) modulation, with amplitude of order H/M - 
and, as a consequence, the amplitude of the modulation 
decreases with increasing k. A plot for the specific case 
of power-law inflation is to be found in and in Fig. 

HI 



the Markov Chain Monte Carlo (MCMC) methods dis- 
cussed below, two widely used approaches have been grid 
searches in parameter space and Fisher matrices, which 
provide theoretical estimates of the likely accuracy that 
can be obtained with a given measurement. A key ob- 
servation that applies to all methods is that we can- 
not simply estimate trans-Planckian parameters after the 
other cosmological parameters are determined. If we add 
new parameters (from trans-Planckian physics or else- 
where) we must re-estimate all the other parameters in 
our model. Failing to do so will result in unphysically 
small error ellipses. 

Grid methods have no intrinsic theoretical shortcom- 
ings, but with more than three or four variables they 
become far too time consuming for a thorough search, 
since the number of points on the grid increases expo- 
nentially with the number of parameters. Fisher matrix 
techniques scale well with the total number of parame- 
ters, but they cannot cope with a likelihood surface that 
contains multiple peaks, and we will see that this is pre- 
cisely the situation we are confronting here. That said, 
we use both methods to make a rough survey of the pa- 
rameter space before discussing our Markov chain results. 

In all cases, we use a modified version of the code 
CMBFAST 4.5.1 with the high-precision option 43] to cal- 
culate the power spectra. The ^-space and /c-space reso- 
lution of CMBFAST was significantly increased in order to 
resolve the trans-Planckian modulations with sufficient 
accuracy. We do not, however, calculate every -^-mode, 
as suggested in 37J — this would render the computa- 
tional cost of exploring such a high-dimensional param- 
eter space prohibitively high, since the computation of a 
single model takes ~ 6 minutes of computing time at the 
resolution used, roughly an order of magnitude longer 
than needed on the same hardware with an unmodified 
version of CMBFAST. However, unlike 37], we are not test- 
ing modulations with independent amplitudes and fre- 
quencies; in our model, very high frequency oscillations 
either have strongly suppressed amplitudes or enormous 
tensor/scalar ratios, or both. Thus, while we need to 
sample in £-space sufficiently to resolve the oscillations 
we are looking for, sampling every £-mode is not a re- 
quirement for the model at hand. 



A. Likelihood Function for an Ideal Experiment 
measuring T,E,B 



III. QUICK SURVEYS 

There are several ways to estimate parameters, and 
their likely uncertainties, in CMB data. In addition to 



^ See Ref. ^3 ^ discussion of the effect of boundary terms 
in the action on this prescripti on. An interesting alternative 
prescription is discussed in Ref. l43 . 



In order to calculate the best possible constraints ob- 
tainable for the trans-Planckian model above, we wish 
to simulate an ideal noiseless CMB experiment with full 
sky coverage. Assuming CMB multipoles are Gaussian- 
distributed, the likelihood function for such an experi- 
ment is given by 

exp [-iDLc-iD,^' 
£ocTT ^ ; 13 
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Here, D^^ is the data vector of spherical harmonic coef- 
ficients 



^Im^ ^fmi ^fm] 



(14) 



where, for example, the temperature map has been ex- 
panded in spherical harmonics as T(n) = Xl^m'^^n*^™' 
and C is the covariance matrix given by 



C = 



Cf^ Cf ^ 
cf^ 



(15) 



The terms Cf^ and Cf^ are zero by global isotropy. 
In the covariance matrix, Cf^ denote theoretical power 
spectra. Now we define the estimator for the power spec- 
tra from the data as: 



XY 



E 
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(16) 



Since the universe is assumed to be isotropic, the likeli- 
hood function is independent of m, and summing over it, 
one obtains, up to an irrelevant constant: 
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c, 



BB 
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BB 
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(17) 



The likelihood has been normalized with respect to the 
maximum likelihood, where Cf^ = Cf^, so that it be- 
haves like a statistic. 

Given a set of observational data points, the goal 
is then to compute the likelihood function for various 
choices of parameters, i.e. a Bayesian analysis. Since 
we are interested in forecasting the performance of fu- 
ture measurements, the "data" is synthetic rather than 
from an actual observation. Various methods for per- 
forming the likelihood analysis have advantages and dis- 
advantages with respect to computational efficiency and 
coverage of the parameter space. In the next section, we 
discuss the results of a grid-based search of the trans- 
Planckian parameter space. 



B. Grids 

Grid calculations are a simple "brute force" method of 
exploring a given parameter space. Given a set of ob- 
servational data points (for example a CMB multipole 



spectrum), the Bayesian best-fit A^-parameter model is 
found by dividing the TV-dimensional parameter space 
into a grid of evenly spaced points and computing the 
likelihood function of each model relative to the data. 
Error bars are assigned by the rate of dropoff of the like- 
lihood function from the best-fit point, which in a coarse 
grid typically involves interpolation between the likeli- 
hoods actually computed on the grid.. Since the number 
of calculations increases exponentially with the number 
of parameters grid methods are inefficient for explor- 
ing large parameter spaces. 

Since we are interested in forecasting the capabilities 
of future CMB measurements we do not use current ob- 
servational data points, but produce "synthetic" data by 
computing a fiducial model representing the assumed un- 
derlying cosmology and assign a set of error bars to the 
points in the fiducial model representing the accuracy of 
the observation. Our hypothetical baseline experiment 
is cosmic variance limited to ^ = 1500 in all four CMB 
spectra: Cj^ , Cj^ , Cf^, and Cf^. While this is an 
unrealistic assumption, it characterizes the amount of in- 
formation which is intrinsically contained in the CMB 
spectra themselves, that is: What is the best we can pos- 
sibly do? Since grid methods are particularly well suited 
to performing a broad survey of a small parameter space, 
we choose our parameter space to be the trans-Planckian 
amplitude H/M^ and the phase (j) of the oscillations, 
which is in principle arbitrary. Other cosmological pa- 
rameters are held fixed. The fiducial model we assume 
for the grid calculation is approximately the best-fit from 
the WMAP data set, = 0.05, 9.c = 0.25, = 0.7, 
h = 0.72, Tri = 0.12, Us = 0.95, dus/dlnk = -0.02. 
We assume a tensor /scalar ratio of r = 0.1, and trans- 
Planckian parameters H/M = 0.022 and (j) = 0.021. 
(These apparently odd choices for the trans-Planckian 
parameters are chosen so that the fiducial model lies ex- 
actly on a grid point in Fig. [21 ) Figure shows the 
result of a coarse grid covering a large range of param- 
eters, showing the substantial degeneracy in the H/M 
- (j) plane, in particular many "islands" in the likelihood 
function representing good fits between the test and fidu- 
cial models. Figure [21 is a finer grid run over a smaller 
region of the same parameter space, showing more detail 
in the shape of the likelihood function. 

Several conclusions can be drawn from this analy- 
sis. First, the projected error ellipses for our hypothet- 
ical cosmic- variance limited measurement are extremely 
small in the H/M plane, indicating a clear ability to de- 
tect the trans-Planckian oscillations in the power spec- 
trum. However, because of the presence of a large num- 
ber of nearly degenerate "islands" in the likelihood func- 
tion, it will be difficult to pin down the precise values 
of both the trans-Planckian parameters, especially the 
phase. This is consistent with the results obtained in 
previous analyses However, we do note that for 

the fiducial model chosen, none of the best-fit "islands" 
is consistent with zero amplitude, H/M = 0. Therefore, 
although the values of trans-Planckian parameters are 
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FIG. 1: A coarse grid (100 points in the range H/M — 
[0,0.05] by 75 points in the range = [0,37r/2]) covering a 
broad region of the trans-Planckian parameter space. The in- 
ner (colored) contours are drawn at the la, 2cr, and 3a levels, 
and the outer (black) contours are at a Ax^ = 50 relative to 
the best fit to better show the shape of the likelihood function. 
The fiducial model has H/M = 0.022 and = 0.021. 



poorly constrained, the case of zero modulation can be 
ruled out to high significance. Figure 01 shows the shape 
of the likelihood surface relative to a fiducial model with 
H/M = 0.01, and Fig. 21 shows the likelihood relative to 
the null hypothesis, H/M = 0, which gives a measure of 
the size of the upper limit one could place on H/M if no 
signal were detected. (We emphasize that the 3cr upper 
limit of about H/M < 0.0025 in Fig. gldoes not take into 
account variations in other cosmological parameters.) 

Why are these degenerate "islands" present in the pa- 
rameter space? In order to understand this result in a 
semi-analytic sense, we consider the effect of the degen- 
eracy on the primordial power spectra. The motivation 
of the following analysis is that likelihood degeneracies 
arise from frequency "beating" between underlying mod- 
ulations in the data and the model. Parameter combina- 
tions where "beating" occurs over the largest /c-ranges 
give rise to discrete likelihood "islands" or local maxima 
in the {H/M) vs (j) plane. 

Defining the trans-Planckian modification factor mul- 
tiplying the standard Bunch-Davies power spectrum as 



A(fc, A,(/),e) = 1 + A 



sm 



2 1 



1 - e A \t 



where A = {H/M)^ we define the parameter / where 



(18) 



J _ /o.ooi Mk, A, (p, e)A(fco, Ap, (pp, €o)dk 
/o°ooi ^(^0, Ao, 00, eo)A(/co, Ao, 00, eo)dk 



(19) 



Here, we have performed the integration over the relevant 
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FIG. 2: A fine grid (100 points in the range H/M = 
[0.020,0.024] by 100 points in the range = [0,0.4]) show- 
ing detail in a small region of the parameter space from Fig. 

The fiducial model is in the center of the plot, and con- 
tours are drawn at la (inner, red), 2a (middle, blue), and 3a 
(outer, green) levels. 
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FIG. 3: Grid results for a fiducial model with H/M = 0.01 
and (/) = 2.0, showing the broadening of the degeneracy for 
smaller amplitude. (Note that the fiducial model is missed 
due to finite grid size effects.) 



scales for the CMB, k - O.OOlMpc"^ - 0.1Mpc"\ If 
/ = 1, there will be a near-perfect degeneracy between 
the parameter sets (Ac^o^^o) and (A,^, e). Assuming 
that e can be constrained by detecting tensor modes, a 
local likelihood maximum is likely to be present at (A, 0) 
even though the underlying "data" has (Ao,0o)- 
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methods in Sec. IIVI 



Q I 1 1 I I III I I I I \ I I \ I I I \ I I I \ I \ \ \ I L 

0.01 0.02 0.03 0.04 0.05 
H/M 

FIG. 4: Grid results for the nuh hypothesis, H/M = 0, es- 
timating an upper hmit on H/M if no trans- Planckian signal 
is detected. 



In the case corresponding to our fiducial models above, 
Ao = 0, only a simple degeneracy 1/A ex irn — (j) exists. 
In Figures [51 we compute / for several values of non- 
zero A = [0.005,0.01,0.04]. While the results in these 
figures are in the analytically tractable power-spectrum 
space. Figure El shows a comparison with Figure El above, 
which was calculated the full Ci formalism. Since con- 
volution with the CMB transfer functions will tend to 
smear out likelihood degeneracies present in primordial 
power-spectrum space, we expect the likelihood "islands" 
to be less prominent in the full Ci calculation than in the 
semi-analytic calculation. 

These results show that the presence of discrete likeli- 
hood maxima in parameter space gives some extra infor- 
mation beyond merely being a numerical headache and 
a barrier to parameter estimation. Their presence in a 
parameter estimation analysis is an indication of an un- 
derlying modulation in the data, and the appearance of 
these degeneracies on the {H/M) vs (j) plane in this model 
is a function of the value of {H/M) of the data. Of course, 
should such an effect be detected in real data, whether 
or not such a modulation is due to a primordial effect or 
some systematic effect in the data, is a question which 
would need to be investigated thoroughly. 

While the grid method is a good tool for exploring the 
global shape of the likelihood function, it is an oversim- 
plification to vary only two parameters in the fit. Error 
bars derived from such a two-parameter analysis will be 
unrealistically optimistic, since a complete analysis must 
take into account degeneracies between the parameters 
of interest and a large number of other cosmological pa- 
rameters. To investigate the effect of including other cos- 
mological parameters in the fit, we apply a Fisher matrix 
technique in Sec. IIII CI and Monte Carlo Markov chain 



Fisher Matrix 



The Fisher matrix technique is a simple and efficient 
method for estimating the expected measurement errors 
when marginalizing over a large number of parameters, 
for which grid methods would be unacceptably time con- 
suming. Measurement uncertainty in cosmological pa- 
rameters is characterized by the Fisher information ma- 
trix aij. (For a review, see Ref. 0.) Given a set of 
parameters {A^}, the Fisher matrix is given by 



EE^c„v-(6.v.ev-,) 



da 



Yl 



I X,Y 



(20) 



where X, F = T, E, B, C and Cov"^ (CxiCyi) is the in- 
verse of the covariance matrix between the estimators 
Cxi of the power spectra. Calculation of the Fisher ma- 
trix requires assuming a "true" set of parameters and 
numerically evaluating the Cxi^s and their derivatives 
relative to that parameter choice. The covariance matrix 
for the parameters {Ai} is just the inverse of the Fisher 
matrix, {a~^).^, and the expected error in the parameter 

Xi is of order \/(cr^^Oii- The full set of parameters {Xi} 
we allow to vary is: 

1. tensor/scalar ratio r, 

2. spectral index n^, 

3. running dus/dlnk^ 

4. normalization A{k^ = 0.002 Mpc"^), 

5. baryon density (^5, 

6. Hubble constant h = i^o/(100 kmsec"^ Mpc~^), 

7. reionization optical depth, th, 

8. Cosmological constant density I^a, 

9. Trans-Planckian amplitude H/M ^ 

10. Trans-Planckian phase (j). 

The total density is fixed such that the universe is flat, 
1^5 + + = 1- Since one needs only to calculate the 
derivative of the likelihood function with respect to each 
parameter. Fisher matrix calculations scale linearly with 
the number of parameters and are therefore very efficient 
for forecasting errors in fits marginalizing over a large 
number of parameters. This technique, however, gives 
no information about the global shape of the likelihood 
surface. Obviously, since we do not know the values of 
the trans-Planckian parameters in the real universe, only 
the size of the errors, not the location of the central value, 
is significant. 
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FIG. 5: (Left) The value of / plotted on the (A,(/)) plane, where the purple (gray) lines, white areas and black areas represent 
/ = 1, / > 1 and I < 1 regions respectively. (Right) The 3-d representation of /. The underlying "data" has 00 = and 
eo = 0.15/16 = e, and the underlying value of Ao ranges from (top to bottom) 0.005, 0.01, 0.04 respectively. 
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0.02 0.021 0.022 0.023 0.024 
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Ao = 0,02, 4)Q = 0.2 




FIG. 6: / plotted on the (A, 0) plane, showing in detail a 
small region for comparison with Figure 13 The underlying 
"data" has = 0.2, eo = 0.15/16 = e, and Ao = 0.02. 



on the value of r in the underlying fiducial model we 
choose, because for small r, the trans-Planckian mod- 
ulation becomes very long-wavelength and thus mimics 
a change in the spectral index ris or running dus/dlnk 
(Fig. [TJ. (For the case r = 0, the "modulation" loses all 
scale-dependence and simply affects the normalization. 
The value of H/M thus decouples from the other param- 
eters determining the shape of the power spectrum, and 
the degeneracy vanishes in this limit.) 



We therefore plot error ellipses in the - H/M plane 
and the dris/dlnk - H/M plane for various choices for 
r running from r = 0.02 to 0.2. The fiducial model as- 
sumed in all cases has H/M = 0.01 and phase (j) = 0. 
Other parameters are identical to those used for the grid 
calculation in Sec. IIII Bl Figure [HI shows the - H/M 
plane for the case of the cosmic- variance limited measure- 
ment, and Fig. [1171 shows the dug/dhik - H/M plane. In 
this case, the measurement is sufficiently accurate that 
the degeneracy between the power spectrum parameters 
n^, dUs/dlnk and the trans-Planckian amplitude H/M 
is negligible, although the ability to constrain H/M is 
significantly degraded for small values of r. These error 
estimates can be directly compared with those generated 
by a full Markov Chain Monte Carlo analysis with simi- 
lar parameters, (Figs. IT^and lT^ . and are in good agree- 
ment. This situation is very different for the case of 
Planck Surveyor, shown in Figs. [Ul and for the - 
H/M and dug/dhik - H/M planes, respectively. We see 
that in this case, the inclusion of a trans-Planckian mod- 
ulation in the parameter set strongly degrades the ability 
to measure the power spectrum parameters. This degen- 
eracy is lifted if the tensor/scalar ratio is large enough, 
effectively disappearing altogether for a tensor/scalar ra- 
tio of order r ~ 0.1. 



We perform a Fisher matrix analysis for two different 
forecast data sets: first, in order to understand the best- 
case scenario, we consider a measurement for which the 
errors are limited by cosmic variance to £ = 1500, sim- 
ilar to that assumed for the grid method in Sec. IIII Bl 
Second, we consider the more immediately realistic case 
of a measurement with error bars similar to those pro- 
jected for the Planck Surveyor satellite, for which we as- 
sume a two-channel experiment measuring both temper- 
ature and polarization. The 143 GHz channel is assumed 
to have an 8-arcminute beam and a pixel sensitivity of 
cTpix = 5.5 /iK, and the 217 GHz channel is assumed to 
have an 5.5-arcminute beam and a pixel sensitivity of 
cTpix = 11.7 /iK. 

We are especially interested in the degeneracy between 
trans-Planckian modulations and the other parameters 
determining the shape of the primordial power spectrum. 
Us and dUs/dlnk. This degeneracy is highly dependent 



Finally, as expected, the measurement uncertainties 
on the trans-Planckian parameters increase significantly 
when a full set of other cosmological parameters is taken 
into account. The Fisher matrix analysis is in good agree- 
ment with the full Markov Chain Monte Carlo analysis 
discussed in Sec. IIVI indicating that a trans-Planckian 
modulation of order H/M = 0.01 is only detectable at 
the 2(7 level in the most optimistic case of a very large 
tensor /scalar ratio r ^ 0.1. In the next section we dis- 
cuss the MCMC analysis in detail and derive the error el- 
lipses relative to a fiducial model with H/M = 0, testing 
our ability to distinguish a signal from the null hypoth- 
esis. This analysis (which includes a longer lever arm of 
i = 2000), indicates that a cosmic- variance limited mea- 
surement can place a 2a upper limit on trans-Planckian 
fiuctuations of around H/M = 0.004 for the case with a 
large tensor /scalar ratio (r = 0.15), and H/M = 0.03 for 
the case with an undetectably small tensor/scalar ratio. 
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FIG. 7: Modulation of the power spectrum for various choices of r: r = 0.02 (top left), r = 0.05 (top right), r = 0.1, (bottom 
left), and r = 0.2 (bottom right). In all cases the solid (cyan) line is the unmodulated spectrum, and the points represent the 
modulated power spectrum as sampled by CMBFAST. 



IV. MARKOV CHAINS 

A. Introduction 

In this section, we use a Markov Chain Monte Carlo 
(MCMC) technique to evaluate the likelihood function 
of model parameters. This approach, proposed by 
in the context of cosmological parameter estimation, 
has become the standard tool for such analyses [e.g. 
m El El S Ei- MCMC is a method to simulate 
posterior distributions. In particular, we simulate ob- 
servations from the posterior distribution V{a\x)^ of a 
set of parameters a given event x, obtained via Bayes' 
Theorem, 



V{a\x) = 



V{x\a)V{a) 
J V{x\a)V{a)da 



(21) 



where V{x\a) is the likelihood of event x given the model 
parameters a and V{a) is the prior probability density. 
The MCMC generates random draws (i.e. simulations) 
from the posterior distribution that are a "fair" sample of 



the likelihood surface. From this sample, we can estimate 
all of the quantities of interest about the posterior distri- 
bution (mean, variance, confidence levels). The MCMC 
method scales approximately linearly with the number 
of parameters, thus allowing us to perform a likelihood 
analysis in a reasonable amount of time. 

A properly derived and implemented MCMC draws 
from the joint posterior density V{a\x) once it has con- 
verged to the stationary distribution. For our applica- 
tion, a denotes a set of 10 parameters. The set is al- 
most the same as used for the Fisher matrix calculation. 
We have 4 "late- time" parameters: the physical energy 
density in baryons /i^, the total physical energy den- 
sity in matter l^rn^^, the optical depth to reionization 
Tri, the Hubble constant in units of 100 km s~^ Mpc 
- h. There are 6 "primordial parameters": the tensor- 
to-scalar ratio r, the spectral slope ng, the running of 
the spectral slope dus/dlnk^ the power spectrum am- 
plitude A, and the parameters of the trans-Planckian 
model defined above, H/M and (j). Event x will be the 
set of observed C^. The priors on the model are that 
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FIG. 8: 2a error ellipses in the Us - H/M plane, calculated 
by Fisher matrix, assuming a fiducial model with Us = 0.93 
and H/M — 0.01 and a cosmic- variance limited measurement 
to i = 1500. The different ellipses show the effect of changing 
the choice of tensor/scalar ratio r in the fiducial model: larger 
values of r make it easier to distinguish the trans-Planckian 
modulation. The top panel shows ellipses for r = 0.02 (black, 
outer), r = 0.05 (blue, next in), r = 0.1 (green, second small- 
est), and r = 0.2 (red, inner). The bottom panel shows the 
same figure zoomed in to show the three inner ellipses. (Neg- 
ative values for H/M are nonphysical, and are included to 
show the full error ellipse.) 



the universe is flat, so that Qm + = 1, and the 
dark energy is a cosmological constant with equation of 
state p = —p. The definitions used for the unmodified 
primordial power spectra of the scalar and tensor per- 
turbations A^j^j^{k) = A:V(27^2)(|7^k|^) oc k""^-^ and 
^Isoif^) = 2ky{27r^){\h^^,\^ + |/ixk|') k^^ are given 
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FIG. 9: 2(7 error ellipses in the Us - H/M plane, calculated by 
Fisher matrix, assuming a fiducial model with ris = 0.93 and 
H/M = 0.01 and an observation with the projected sensitivity 
of the Planck Surveyor satellite. The different ellipses show 
the effect of changing the choice of tensor/scalar ratio r in the 
fiducial model: larger values of r make it easier to distinguish 
the trans-Planckian modulation. The top panel shows ellipses 
for r = 0.02 (black, outer), r = 0.05 (blue, next in), r = 
0.1 (green, second smallest), and r = 0.2 (red, inner). The 
bottom panel shows the same figure zoomed in to show the 
three inner ellipses. 



in Section ITTI 



B. Methodology 

The experiment in the case study is considered to mea- 
sure TT, TE, EE, BB power spectra to ^ = 2000; beyond 
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FIG. 10: 2(7 error ellipses in the dris/dlnk - H/M plane, 
calculated by Fisher matrix, assuming a fiducial model with 
ris — 0.93 and H/M — 0.01 and a cosmic- variance limited 
observation. The different ellipses show the effect of changing 
the choice of tensor/scalar ratio r in the fiducial model: larger 
values of r make it easier to distinguish the trans-Planckian 
modulation. The top panel shows ellipses for r = 0.02 (black, 
outer), r = 0.05 (blue, next in), and r = 0.1 (green, smallest). 
The bottom panel shows the same figure zoomed in to show 
the three inner ellipses. 



this, the systematic errors in removing the contamination 
from secondary CMB anisotropics are likely to be non- 
negligible. We include the effect of weak lensing on all 
power spectra, since this is the "theoretical limit" in the 
accuracy of the measurement - weak lensing aliases power 
between multipoles and tends to smooth out small-scale 
oscillations in the data. Therefore, any of the further 
uncertainties that must be accounted for in a practical 
CMB experiment, such as foreground removal, partial 
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FIG. 11: 2a error ellipses in the dus/dlnk - H/M plane, 
calculated by Fisher matrix, assuming a fiducial model with 
Us = 0.93 and H/M = 0.01 and an observation with the 
projected sensitivity of the Planck Surveyor satellite. The 
different ellipses show the effect of changing the choice of ten- 
sor/scalar ratio r in the fiducial model: larger values of r make 
it easier to distinguish the trans-Planckian modulation. The 
top panel shows ellipses for r = 0.02 (black, outer), r = 0.05 
(blue, next in), r = 0.1 (green, second smallest), and r = 0.2 
(red, inner). The bottom panel shows the same figure zoomed 
in to show the three inner ellipses. 



sky coverage etc, will degrade the constraints obtainable 
on the trans-Planckian parameters relative to those found 
in this study. 

For the simulated "data", we consider two cases: one 
with significant primordial tensor modes (denoted Case 
H), and one with an undetectably small contribution of 
primordial tensors (denoted case L). Case H is a mass- 
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term inflation model with tensor-to-scalar ratio r = 0.15 
which is consistent with WMAP data at the la level. 
Case L is simply the WMAP LCDM concordance model 
with r = 0.00013 (quadrupole ratio /C^ - 6 x 10"^). 

The objective of this study is to see how far from the 
null case the trans-Planckian parameters must be before 
they are distinguished from the null model by at least 2a. 
Consequently, these flducial models do not contain any 
trans-Planckian oscillations. Our flducial models also do 
not contain any running of the primordial scalar spectral 
index; however, the effect of this parameter would be 
degenerate with a very long wavelength trans-Planckian 
modulation. Hence we marginalize over this parameter 
in our analysis. 

In the practical implementation of the MCMC, we ran 
16 chains for Case H and 32 chains for Case L, each 
started at random 3a steps away from the flducial model. 
By making use of preliminary test chains, the proposal 
distribution was calibrated using the covariance between 
the parameters in order to minimize correlations; the 
chains were conflgured to step in the principal directions 
of the parameter covariance matrix optimal ~ 24% accep- 
tance rates for case H; for case L, only a few % acceptance 
rates were achieved even after several recalibrations of the 
proposal matrix, indicating the complexity of the likeli- 
hood surface. The chains were run on a high-performance 
64-bit Opteron Beowulf cluster, and 300,000 and 500,000 
total models were computed for Case H and Case L, re- 
spectively. The presence of "islands" in the likelihood 
space meant that the chains converged relatively slowly, 
compared to what is typical for a flt to CMB observations 
when the trans-Planckian terms are absent. In fact, we 
stopped the chains for Case L slightly short of full con- 
vergence, but long after the predicted parameter regions 
had stabilized. 



C. The Shape of the Likelihood Surface 

Tables HI and ITTI show the flducial models and recovered 
marginalized 1-d 68% probability parameter constraints 
on the ten model parameters for Case H and Case L, 
respectively. Figures and show the marginalized 
2-d joint 68% and 95% confldence regions for pairs of 
the six "primordial" parameters for the same models. In 
Case H, which has potentially observable primordial ten- 
sor modes, a deviation of H/M from null is detectable at 
greater than 2a level if {H/M) > 0.004 in an ideal exper- 
iment, which renders H/M ~ 0.01 potentially detectable 
in a real- world experiment. In Case L, which has un- 
observably small primordial tensor modes, one requires 
{H/M) > 0.03 in order to be detectable at greater than 
the 2a level in an ideal experiment, rendering such a de- 
tection implausible in the real world. Ref. discussed 
a degeneracy between e and {H/M) which appears to 
be essentially eliminated in our model due to the extra 
constraint on e from tensor modes. 

The key observations to be drawn from these results 



TABLE I: Case H: A fiducial model with r = 0.15 and no trans- 
Planckian modulations. The input model, recovered mean param- 
eters and the 68% probability level of the 1-d marginalized likeli- 
hood. The derived parameter erg (the linear theory prediction for 
the amplitude of fiuctuations within 8 Mpc/h spheres) is given in 
the last row. 



Parameter 


Input 


Recovered 


O 

\Lbtl 


U.Uzz / o 


U.Uzz /o ± U.UUUUd 




0.3097 


0.3096 ± 0.0016 


h 


0.684 


0.684 ± 0.001 




0.112 


0.112 ±0.002 


ns{k^ = 0.002 Mpc-^) 


0.962 


0.964 ± 0.013 


dus /d\nk 





0.000 ± 0.004 


A{k^ = 0.002 Mpc-^) 


0.879 


0.878 ± 0.016 


r{k^ = 0.002 Mpc-^) 


0.150 


0.150 ±0.004 


H/M 





0.002 ± 0.002 







(degeneracy; mean 0.02) 






0.871 ±0.001 



TABLE II: Case L: A fiducial model with r = 0.00013 and no 
trans-Planckian modulations. The input model, recovered mean 
parameters and the 68% probability level of the 1-d marginalized 
likelihood. The derived parameter erg (the linear theory prediction 
for the amplitude of fiuctuations within 8 Mpc/h spheres) is given 
in the last row. 



Parameter 



Input 



Recovered 



Qbh^ 0.02380 

Qm 0.2943 

h 0.700 

Tri 0.166 
ns{k^ = 0.002 Mpc-^) 0.994 

dus /d\nk 

A{k^ = 0.002 Mpc-^) 0.900 



0.02379 ± 00006 
0.2945 ± 0.0016 
0.700 ± 0.001 
0.166 ±0.002 
0.993 ±0.010 
0.000 ± 0.003 
0.900 ± 0.020 



r{k^ = 0.002 Mpc-^) 0.00013 0.00014 ± 0.00004 

H/M 0.02 ± 0.02 

(degeneracy; mean 0.05) 
as 0.927 ± 0.001 



are: 



1. The constraint on H/M differs by an order of mag- 
nitude between Case H and Case L. This reinforces 
the point made above that if our Universe happens 
to be in a favorable part of the parameter space 
where the energy scale of inflation is high enough 
to produce an observable amplitude of primordial 
tensor modes, then one has a significantly improved 
chance of detecting any trans-Planckian modifica- 
tions to the spectrum whose properties are similar 
to those of the case tested here. 

2. The cosmological parameters are constrained at the 
~ 1% level except for the scalar spectral index and 
the amplitude of fluctuations at /c = 0.002 Mpc~^, 
which are only constrained at the ~ 10% level. 
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FIG. 12: Case H: A fiducial model with r — 0.15 and no trans-Planckian modulations. The 2-d joint 68% (dark) and 95% 
(light) confidence regions for pairs of parameters in the set {r, ris, dris/dhik^ A, H/M, 0} obtained via MCMC methods. In 
the case of each parameter pair, all the other eight parameters in the model have been marginalized over. The discrete "blobs" 
in the (ris, dus/ dink) plane are not physical but a graphical effect arising as consequence of binning of a narrow likelihood 
degeneracy. 



Therefore, another cosmological data-set that aids 
in constraining these parameters (for example, the 
proposed SKA survey, 50]) will minimize the de- 
generacies with these parameters shown in FigsIT^ 
and El and improve parameter constraints further 



V. CONCLUSION 

We believe that this analysis is the most thorough in- 
vestigation of the detectability of a trans-Planckian mod- 
ulation to the primordial power spectrum that has been 
performed to date. We consider three complementary 
approaches - a simple grid search, a Fisher matrix eval- 
uation of the likely error ellipses, and a Monte Carlo 
Markov Chain fit to a simulated CMB spectrum, and 



find that they are all in broad agreement. Moreover, we 
also explain the "islands" seen in the likelihood space 
in previous papers on this topic |3 US ? show that 
their distribution and properties can be understood and 
reproduced via a simple analytic argument. 

The approach we have taken here is analyze a specific 
ansatz for introducing a minimum length into the calcu- 
lation of the perturbation spectrum in a general model of 
slow-roll inflation 11]. We do not claim that this model 
is a correct description of the trans-Planckian contribu- 
tion to the perturbation spectrum, but intend it as a 
case study of what might be possible if one has a spe- 
cific and well-motivated correction to the spectrum that 
leads to a /c dependent modulation. There is consider- 
able theoretical uncertainty surrounding this point, and 
it is entirely possible that a future rigorous calculation of 
the perturbation spectrum within string theory or some 
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FIG. 13: Case L: A fiducial model with r — 0.00013 and no trans-Planckian modulations. The 2-d joint 68% (dark) and 95% 
(light) confidence regions for pairs of parameters in the set {r, ris, dris/dhik^ A, H/M, 0} obtained via MCMC methods. In 
the case of each parameter pair, all the other eight parameters in the model have been marginalized over. The discrete "blobs" 
in the (ris, dus/ dhik) plane are not physical but a graphical effect arising as consequence of binning of a narrow likelihood 
degeneracy. 



other model of ultra high energy physics will predict that 
there is no modulation to be observed, and the bound 
calculated here will not apply. 

We believe that the constraint obtained here is appli- 
cable to a general class of modulated spectra. In practice, 
one can imagine adding an arbitrary modulation to the 
primordial spectrum, as considered by Okamoto and Lim 
J35] , in which case one has parameters which describe the 
amplitude, wavelength, and phase of the modulation as 
a function oik. In our model there are only two free pa- 
rameters - the phase of the modulation, and H/M which 
determines the amplitude. The phase will always be ar- 
bitrary in the absence of detailed theory that matches 
specific scales in the present universe to those during in- 
flation, which would require a full understanding of the 
post-inflationary expansion history and the physics of re- 
heating, while the amplitude is directly related to the 



scale of new physics, which is what we are trying to mea- 
sure. However, the wavelength of the oscillation turns 
out to be given in terms of the slow-roll parameter e, 
which is itself directly related to the tensor /scalar ratio 
r. This dependence is easy to understand, since the form 
of the modulation is not determined by the (fixed) bound- 
ary condition, but by the slow variation in the horizon 
size H~^ with time. The precise functional form of the 
connection between e and the modulated spectrum may 
vary in other specific models. However, it is very likely 
that the modulation can always be expressed in terms of 
H/M and the scale dependence of - or thus e, and 
perhaps the higher order slow-roll parameters, which are 
also reachable through their contribution to the scalar 
spectral index. Consequently, we believe that the re- 
sults we have seen here will generalize to other modu- 
lated spectra, even though we have reduced the number 
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of trans-Planckian parameters relative to those consid- 
ered by Okamoto and Lim. 

The second general conclusion we draw is that the de- 
tectability of any modulated spectrum depends on the 
value of H. Since the wavelength of the modulation de- 
pends on e, detecting the primordial tensors provides an 
orthogonal constraint on the value of this parameter, and 
thus constrains the trans-Planckian corrections. Since H 
fixes the energy scale at which inflation occurs, it is ar- 
guably the single most interesting cosmological parame- 
ter that is not currently fixed by observations. The analy- 
sis here simply adds to the importance of this parameter, 
since we have shown that measuring H will put tighter 
bounds on any trans-Planckian corrections to the spec- 
trum. 

Since the theoretical uncertainty about the form and 
existence of these corrections is yet to be resolved, we 
cannot advocate mounting an observational campaign 
solely to look for this kind of signal. However, there is al- 
ready good and sufficient reason to make an all out effort 
to measure H during inflation, via high precision mea- 
surements of the CMB and particularly its polarization. 
We believe that mission planners will want to be aware 
of the possibility that precision CMB measurements may 
potentially probe physics at the string scale, since this 
will not add to the cost of any mission but represents 
an exciting and additional use of the data they can be 
expected to return. 

In quantitative terms, we see that the with a large 
(r ~ 0.15) value of the tensor to scalar ratio, a "per- 
fect" map of the primordial temperature and polariza- 
tion anisotropics to the CMB could rule out a trans- 
Planckian modulation to the spectrum at the 2cr level 
for H/M ~ 0.004. This is roughly the same as the result 
reported by Okamoto and Lim. We vary a larger set of 
cosmological parameters than they do, which reduces the 
level of precision we can hope for. However by adding the 
tensor contribution to the analysis and removing a free 
parameter from the specification of modulation we can 
tighten the bound. Like Okamoto and Lim, we consider 
a perfect measurement of the CMB and any real experi- 
ment will be contaminated by imperfectly removed fore- 
grounds. Consequently, we believe that a value of H/M 
on the order of 0.01 represents a reasonable lower bound 
on what can be detected in practice, provided the am- 
plitude of the tensor contribution is significant. On the 



other hand, adding other orthogonal datasets, such as a 
measurement of h and the late time expansion history of 
the universe from a SNAP-like mission, or constraints on 
the primordial spectrum from future large scale structure 
surveys will further constrain the error ellipses relative to 
those which can be achieved via the CMB alone. 

At this point, the theoretical uncertainty surrounding 
this calculation makes it unwise to push this analysis sig- 
nificantly further. However, we believe that consider- 
able progress has been made in understanding the physi- 
cal issues surrounding the trans-Planckian corrections to 
the primordial spectrum over the last few years. Conse- 
quently, it is not excessively optimistic to hope that the 
remaining issues can be resolved on a shorter timescale 
than it will take to perform a measurement of the CMB 
that even approaches the precision we have assumed in 
our analysis. In this light, we are particularly interested 
in the effective field theory approach taken in T^, and 
extended in the more recent preprint 52]. In this case, 
the operators which can contribute to a modification to 
the spectrum can be cataloged within effective field the- 
ory. From the observational perspective, an analogous 
calculation to the one presented here will constrain the 
values of the prefactors in front of these operators, and 
we intend to pursue this in future work. 
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